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Abstract 

To constrain cosmological parameters, one often makes a joint analysis with different combina- 
tions of observational data sets. In this paper we take the figure of merit (FoM) for Dark Energy 
Task Force fiducial model (CPL model) to estimate goodness of different combinations of data 
sets, which include 11 widely- used observational data sets (Type la Supernovae, Observational 
Hubble Parameter, Baryon Acoustic Oscillation, Cosmic Microwave Background, X-ray Cluster 
Baryon Mass Fraction, and Gamma- Ray Bursts). We analyze different combinations and make a 
comparison for two types of combination based on two types of basic combinations, which are often 
adopted in the literatures. We find two sets of combinations, which have strong ability to constrain 
the dark energy parameters, one has the largest FoM, the other contains less observational data 
with a relative large FoM and a simple fitting procedure. 
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I. INTRODUCTION 



Since the discovery of cosmic acceleration expansion of the universe l|, |2| , its origin has 
been a hot topic in modern cosmology and theoretical physics. Various dark energy (DE) 
models explaining the acceleration expansion have been proposed 3l4l4j|. which originate 
from extra dimensions, string theory, supergravity, extended standard model of particle 
physics, and so on. Most of the DE models can fit astronomical observational data. At 
present, the main task of DE studies is to constrain the present equation of state (EoS) 
of DE Wdeo, the fractional energy density Qdeo and to give the possible evolution form of 
Wde{z). Several parameterization forms of WdJz) have been proposed, such as the CPL 



parameterization Wde{,z) = Wq + Wiz/{1 + z) [l5|, |l6| and the parameterization Wde{,z) 



Wo + WiZ 



17H20|. Fitting the parameterization models of DE with present observational 



data, constraints of DE can be obtained. The fitting results favor that the present value of 
Wde is very close to —1 and its variation is very slow. But it still not be able to confirm 
whether Wde is evolving or just a constant. Of course, these constraints are also dependent 
on the parameterization forms used. As a result several methods have been proposed to get 
model-independent constraints of DE from observations [21 



27l |. such as the uncorrelated 



bandpower estimate (UBE) 28| and the spline method 291431 1. 

On the other hand, the fitting results also depend on the observational data used and there 
might exist "tensions" among the estimated EoS parameters from different data sets 
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A single observation still gives poor constraints on the parameters of cosmological models. 
For example. Type la supernovae observation, which has large number of datapoints, gives 
poor constraints on Vtmo and Wde because of the degeneracy between fi^o and Wde[39|. To 
break the degeneracy and get good constraints of the parameters, one always combines dif- 
ferent types of observation to fit cosmological models. It is not surprising that fitting a 
DE model with different combinations of observational data sets, one would get different 
constraints on DE. There seemingly does not exist a general rule to combine different obser- 
vations for fitting the cosmological parameters. One uses different combinations in different 
works. 

In this paper we would like to make a comparison for different combinations of widely- 
used data sets for DE, which contain data of Type la Supernovae (SN), Observational 
Hubble Parameter (HUB), Baryon Acoustic Oscillation (BAO), Cosmic Microwave Back- 
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ground (CMB), X-ray Cluster Baryon Mass Fraction (CBMF), Gamma- Ray Bursts (GRB), 
etc. The fitting properties of observational data will be studied by using the Dark Energy 



Task Force (DETF) |40l442| fiducial model (GPL model), i.e., all combinations of data will 
be fitted with the GPL parameterization. The best-fitted values of ^^0,^0,101 and their 
68%, 95% G.L. errors for each combination will be given. To estimate constraints on Wde and 
compare the goodness of data, we will calculate the figure of merit (FoM) for each 

combination, which is proportional to the inverse area of the error ellipse in the wq ~ wi 
plane. Namely we will take the FoM as the diagnostic to quantify the goodness of each 
combination of observational data to constrain the DE parameters. 

The paper is organized as follows. In Section II, we introduce 11 sets of observational 



data and study their fitting properties on the basis of Type la supernovae Union2 data |43|. 
In Section III we study the effect of different combinations of data sets on constraining DE 
parameters. We divide the combinations into two types, based on two basic combinations, 
SN+BAOI+GMB and SN+A+R, respectively. The properties of various combinations will 
be analyzed and compared, and the best combination (which has the largest FoM) will be 
given. Main results are summarized in Section IV. 

II. OBSERVATIONAL DATA SETS 

One of important astronomical observations to constrain DE is the distance measurement 
(such as SN and shift parameter R). To constrain DE models with these data, one must give 
the form of Hubble function H{z) (or E{z) = H[z) / Hq). In a flat Friedmann-Robertson- 
Walker (FRW) universe, one has 

E\z) = ^],o(l + Zf + fibo(l + Zf + ^]rf„o(l + zf + ^deoF{z), (1) 

where fij-o, ^dmo and VLdeo are present values of the fractional energy density for radia- 
tions, baryons, dark matter and dark energy respectively, with Vtro + + ^dmo + ^deo = 1- 
The energy densities of baryons and dark matter are often written together as VL^q + fidmo = 
VLmo- The radiation density is the sum of photons and relativistic neutrinos |44|: 

=fif (l + 0.2271Ar„), (2) 

where is the number of neutrino species and = 2.469 x lO'^/i-^ for T^^^, = 2.72bK 
{h = Hq/100 Mpc ■ km ■ s~^). The evolving function F{z) for DE depends on the equation 



of state, Wde{^)j of DE models, 

= e^^ (3) 



We will fit observational data with the CPL parameterization. In that case, one has 

F{z) = (1 + ^)3(i+-o+-i) exp[-|!^^]. (4) 

When fitting the CPL model with data of distance measurement, the parameters fimo (or 
^deo))'W^O) "W^i; and h are needed. While for some other data, such as those of CMB and BAO, 
one needs to calculate physical quantities other than the distance (such as the redshift of 
decoupling, Zg, and the comoving sound horizon at decoupling, r^) and the parameter f2f,o 
must be involved. Whenever the fitting procedure needs to use fifeo? oii^ more datapoint will 
always be adopted to constrain f2f,o and h |45|: 

Qfeo^^ = 0.022 ± 0.002 . (5) 

In addition, in all calculations we will assume the prior that the age of the universe Tq 
satisfies lOGyr < Tq < 20G?/r. 

We will use Union2 SN data 43| and other widely-used 10 observational data sets to 



constrain the CPL model. When studying the fitting properties of these 10 observational 
data, we will combine each data set with the SN data to fit CPL model. Since we are only 
interested in the constraints on DE, only the best-fitted values of VLmo-,Wo, Wi and their 68%, 
95% C.L. errors will be given, which are obtained by using the Markov Chain Monte Carlo 
method. To estimate and compare the goodness of constraints from the observational data 



sets, we calculate the figure of merit (FoM) 40l-l42l. |46| for each combination of data sets. 



which is proportional to the inverse area of the error ellipse in the wq ~ wi plane: 

FoM = [detC(wo,wi)]"^/^ , (6) 

where C{wq,Wi) is the covariance matrix of Wq and Wi after marginalizing out all of the 
other cosmological parameters. Larger FoM means stronger constraint on the parameters 
since it corresponds to a smaller error ellipse. 



A. Type la Supernovae (SN) 



Type la supernovae observation gives the direct evidence for the existence of the acceler- 
ation expansion of the universe. In this work we take the Union2 SN dataset 43|, for which 
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data 




Wo 






FoM 


SN 


n /1 1Q+0. 090+0.133 


n QR+O.46+1.19 


c: CI +6.99+7.59 
'-'•c^ 7Q— 9zL HQ 


541.431 


0.647 


SN+GRB 


n /I iq+0. 091+0.132 
u . uoo u.o ±u 


n OO+0.42+1.16 

U.OU U.O'-t 


r n9+6.48+7.13 
'-'•^■^ — 7 '^9— 1Q 


564.727 


0.700 


SN+CBMF 


n 070+0.085+0.126 

U.O / On noR—O 097 


r, Q7+0.31+0.71 
"•>^ ' —0 23-0 47 


9 77+2.93+4.36 
^- ' 1 —4 qn— 1 3 nq 


583.725 


1.356 


SN+HUB 


n 981 +0.134+0.170 
U.^01_Q Q27_o.217 


1 nn+0-22+0.47 
^•'^^-0.12-0.34 


n 91+1.54+1.98 
^•^^-2.82-8.46 


550.391 


2.315 


SN+CMB 


n 97/1+0.025+0.051 
^■^'^-0.017-0.036 


1 ni +0-16+0.29 
^■'^-'^-0.14-0.34 


n ni +0.73+1.45 
"-'•"-'^-0.84-1.54 


542.697 


21.688 


SN+R 


n 977+0.021+0.045 
'-'■^' ' -0.021-0.040 


1 ni +0-17+0.31 
-'-■^-'--0.14-0.30 


n nfi+0-75+1-38 
^•^"-0.84-1. 64 


542.642 


21.376 



TABLE I: The best-fitted values with 68% and 95% C.L. errors of J^mO; ^O; ^ii and figure of merit 
for different combinations of data sets in the CPL model. 



xIn is given by 
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2 _ \^ [^J'th{z^) — l^ob{Zi)]'^ 



4 = 1 



(7) 



where the theoretical distance modulus ^ith is defined as 



l^th{z) = 5 logio Dl + fJ^o , 1^0 = 42.384 - 5 log^o h , 

where the luminosity distance = (l + z) dx/E{x). One can expand Eq.(I7]) with respect 
to jj,Q as 

xIn = ^ + + (8) 



where 



557 



i=l 
557 



5 = E 



[fith{zi; fiQ = 0) - fiob{zi)f 
o\z,) 

jJ^thizf, jJ^O = 0) - fiob{Zi) 



i=l 
557 



cr^iz,) 



(9) 



^ = e4 



i=l 



cr^iZi 



The XsN a minimum as 

xIn = A- B'/C . (10) 

We will use xIat instead of Xsn^ ^^is way the nuisance parameter /ig is reduced. It is 
equivalent to performing an uniform marginalization over /ig. 

In Table [U it shows that the SN data give a largish ^.^o and the value of FoM is rather 
small, which means that the constraint on the parameters is weak. This is mainly due to the 



degeneracy between and Wde- Thus the SN data are often combined with other types 
of observations, especially that can constrain Q^o tightly, such as BAO and CMB. 



B. Gamma-Ray Bursts (GRB) 

Gamma-Ray Bursts (GRB) are the most intense explosions we have observed in the 
universe, most GRB observed are in the range 0.1 < z < 8.1. Thus, GRB are considered to 
be a complementary probe to SN and a hopeful probe at high redshift. Here we adopt the 



59 Hymnium GRB at 2; > 1.4 from 



he 109 GRB 



48|. 



which are obtained by calibrating t 
with Amati relation, using the cosmological independent method proposed in 

The operation of GRB datapoints is the same as that of SN. To fit those data with DE 
models it should be combined with SN data, i.e., Xgrb should be calculated together with 



2 

XsN 



where 



x'sN+aRB = A- B'lC . (11) 



[fithiZj] /io = 0) - fiob{Zj)f 



'3) 



557+59 f „x f . 



j=l ^ 1' 

557+59 



The datapoints in the summation contain 557 SN and 59 GRB. 

The results for SN+GRB are shown in Table I, which are almost the same as that of SN. 
With addition of 59 GRB data, the increase of (about 23) is much less than the increase 
of number of datapoints, which implies that the present precision of the GRB observation 
is not small enough. As a result it means that adding GRB data to SN data is not quite 
helpful to constrain the DE model. 

C. X-ray Cluster Baryon Mass Fraction (CBMF) 

The baryon mass fraction in clusters of galaxies (CBMF) are also used to constrain 



cosmological parameters. Here we adopt the updated 42 observational f„as data in 



49|. 



The X-ray gas mass fraction is defined as the ratio of the X-ray gas mass to the total mass 
of a cluster, which should be a redshift-independent constant. With the ACDM reference 
cosmological model, one has 

( m3/2 



jjACDMf 

Da(z) 



(13) 



1 + S 

where A is the angular correction factor, which is due to the change in angle subtended by 
^^2500 as the underlying cosmological model is changed: 

H(z)Da(z) 



(14) 



V ^2500 ; \[H{z)Da{z)]^^'''', 

Here we adopt the best-fitted average value t] = 0.214, and the proper angular diameter 
distance is given by 



Dl{z) 
{1 + zy 



(15) 



For other parameters K,7, b, s in Eq.( JT3l) . please refer to Refs. 49|, |50|. Here we introduce 
a new parameter A 51| : 



1 s fi. 



then one has 



Next one can expand 



fgas{Zi) = XA 



D 



AC DM I 



DAiZi) 



3/2 



— Xfgas(^Z{ 



2 _ {fgas{Zi) fgas,ii 

XCBMF = 



with respect to A and obtain 



where 



XcBMF = AX^ - 2BX + C, 



^ _ ^ fgas{Zi) 



(16) 
(17) 
(18) 

(19) 



a. 



2 ' 



^ _ fgas{Zi)fgas,i 



(20) 



f fobs \2 
\J gas,i) 



a? 



As in the case of SN, instead of XcBMFy here we adopt the minimum of Xcbaif with respect 
to A, 

x'cBMF = C- By A , (21) 



as for CBMF. 

As shown in Table I, the constraint from SN+CBMF is also weak, but better than that 
of SN+GRB. The best fitted value of VL^q is still largish. The addition of 42 CBMF to SN 
data does not improve much. 



D. Observational Hubble Data (HUB) 



The observational Hubble data can be obtained by using the differential ages of passively- 
evolving galaxies as 

H^--^^. (22) 



We use 12 observational Hubble data from 
square is defined as: 



z At 



52j and 53|], which is tabled in 50|]. The chi- 



XhUB - 2 
i=l » 



(23) 



In this case we use these values of H{zi) directly, which can break the degeneracy between 
Qmo and Wde- With these data, as shown in Table I, there is an increase of FoM and the 
best-fitted ^2^0 becomes much smaller than in the case of the SN data used only. 



E. Cosmic Microwave Background (CMB) 

1. CMB 

In the CMB measurement, the distance to the decoupling epoch can be accurately de- 
termined from the locations of peaks and troughs of acoustic oscillations. For simplicity, 
we use three parameters la, R, Zg obtained from WMAP7 (rather than the full data of 



WMAP) 5J] to constrain cosmological models. Here Zg is the redshift of decoupling 55| . 



the shift parameter R is the scaled distance to the decoupling epoch: 

and la is the angular scale of the sound horizon at the decoupling epoch: 

rjas) 

la = , (25j 

rs{as) 
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where r{z) = dx/H{x) is the comoving distance and rs(as) is the comoving sound horizon 
at the decouphng epoch: 



a'^H{a) 



da, a. 



(26) 



where the sound speed Cs{a) = 1/^3(1 + Rba) and Rh = 30^0/4^2^0 is the photon-baryon 
energy density ratio. 

The of the CMB data is constructed as: 



where 



XCMB — ^'^Cjy^X 
( 1 _ Qfio no \ 



(27) 



X 



\ 



la- 


302.09 


R- 


- 1.725 


z. - 


1091.3 



(28) 



and the inverse covariance matrix 



^-1 



V 



2.305 29.698 
29.698 6825.270 
-1.333 -113.180 



/ 

-1.333 
-113.180 
3.414 



(29) 



In this case the parameter VL^q is involved. From the definition of the shift parameter R 
[see (p4|) ]. it is obvious that fimo can be well constrained, and the degeneracy between VLmo 
and Wde could be broken in some sense. As shown in Table HI the FoM and constraints of 
parameters have a great improvement from those in the case of SN data used only. Therefore 
the CMB data are a good supplementary of SN data in constraining DE models. 
2. Shift Parameter R 

As an alternative, one can only adopt the shift parameter R from WMAP7 Q] to con- 
strain the cosmological parameters. In this case the corresponding is defined as 



i?- 1.725 
0.018 



(30) 



It can be seen from Table [T] that the fitting results of SN+R are very close to those of 
SN+R+/a + (i.e., SN+CMB). Note that in this case one needs not to handle with fi^o and 
the calculation is much simpler. The addition of a single datapoint R to SN data greatly 
improves the fitting results, because R can alleviate the degeneracy between Vtrno and Wde 
well. 
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data 








y2 • 


FoM 


SN+BAOI 


n /ion+0-036+0.068 

U.4i:Zu_r\ Ofifi— 177 


n O/I+0.23+0.71 
u.o^ — u.oo 


c- 7Q+6.65+7.06 
o. 1 y_o fio— in 9^ 


542.186 


1.471 


SN+BAOII 


n /19O+0.083+0.122 
" " — 033—0 359 


n sn+0-46+1.20 

^•O^— 35-0 59 


fi ^n+7-77+8.45 
u.<JU_o 7a_26 82 


542.643 


0.632 


SN+BAOIII 


n /1 91 +0.035+0.066 
U.*^i_0.066-0.190 


n 89+0.23+0.67 
^•0^-0.33-0.55 


cr 8Q+6.73+7.13 
0.oy_2 i4_8 69 


542.134 


1.514 


SN+A 


n 978+0.025+0.047 
"-'•^'0-0.018-0.037 


1 ni +0-17+0.33 
-'-•^-'--0. 13-0. 28 


n in+1-04+1.89 
U.lU_i 27-2.57 


542.643 


13.223 


SN+RBAO 


n oQ/i +0.026+0.056 
U.oy^_Q Qg]^_Q 241 


n Qo+0. 18+0.48 
U.yo_Q 2i_o.41 


q 78+4.75+5.05 
"-••'0-0.038-4.35 


541.602 


2.586 



TABLE II: The best-fitted values with 68% and 95% C.L. errors of ^mOi wq, wi, and figure of 
merit for different combinations of data sets in the CPL model. 

F. Baryon Acoustic Oscillation (BAO) 



As baryons and photons are tightly coupled at early times, gravity and pressure gradients 
induce an acoustic oscillation in the baryon-photon fluid. Since the baryonic matter interacts 
gravitationally with the dark matter, the acoustic oscillations leaves some fingerprint in 
the matter power spectrum. The BAO peak length scale is set by the sound horizon at 
decoupling, ~ 10^ Mpc. We now introduce several types of BAO data. 
1. BAOI 

The first one is the BAO distance measurements obtained at z = 0.2 and z = 0.35 from 



joint analysis of the 2dFGRS and SDSS DR7 data 



Dv{0.2) 



0.1905 ±0.0061 , 
0.1097 ±0.0036 , 



DviO.35) 

where Tsizd) is the comoving sound horizon at the baryon drag epoch Zd 



and 



Dviz) 



dx 



1/3 



(31) 
(32) 

(33) 



,0 H{x)} Hiz) 

encodes the visual distortion of a spherical object due to the non Euclidianity of a FRW 
spacetime. The Xbaoi given by 



xIaoi = X^V-'X 



where 



X 



Dv{0.2) 

rs(zd) 
DviO.35) 

10 



-0.1905 
-0.1097 



(34) 



(35) 



and the inverse covariance matrix 

( 30124.1 -17226.9 , , ^ 

V-'=\ \ ■ (36) 

\^ -17226.9 86976.6 

With addition of the BAOI data to SN data, there is a bit improvement of FoM, as shown 
in Table [Tll The best-fitted parameters are also very close to those from the case of SN data 
used only. This implies that BAOI data cannot alleviate the degeneracy between Vtrno and 

Wde- 

2. BAOII 

Using two distance measurements of BAOI, one can derive a model independent BAO 
distance ratio 56| 

^.1.736^0.085. (37) 

which is independent of VLi,q. For this single datapoint, one defines 

_ / Dy(0.35)/Z}y(0.2)- 1.736 V 
Xbaoii-^ 0065 ) ■ ^^^^ 

The FoM of SN+BAOII is even smaller than that of SN data, and the Xmm larger than 
that of SN+BAOI data. 

3. BAOIII 

In this case we combine the data of BAOII with the BAO distance measurements at 
z = 0.275 



^'^^'^'^ - 0.1390 ±0.0037 . (39) 



^^(0.275) 

These two datapoints have been widely-used in the literatures. In this case the is defined 

as 

Xbaoiii = 7^7^ + 7777^ • (^O) 



As shown in Table II, the results for the combination SN+BAOIII are almost the same as 
those of SN+BAOI. 

4. Distance Parameter A 

The distance parameter A is often used in the literatures, which is the measurement of 
BAO peak in the distribution of SDSS luminous red galaxies jssj l: 



A = n^/2^(0.35)-i 



/3 





r^-^^ dz ' 


0.35 J 


E{z)\ 



2/3 



(41) 
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This quantity is independent of Qbo- The value of A is determined to be 0.469(?2s/0.98) "^-^S-i- 
0.017, where Ug = 0.963 is the scalar spectral index, which has been updated from the 
WMAP7 data, is defined as 

_ {A - 0.472)^ 
Xa- qq^72 • (42) 

The supplementary of the single datapoint A to SN data gives a large increase of FoM, 
because Qmo can be well constrained and the degeneracy between Q^o and Wde is alleviated. 
5. RBAO 

SDSS data can also be used to measure the radial (line-of-sight) baryon acoustic scale, 
which is independent from the previous BAO measurements which were averaged over all 
directions or in the transverse direction. In that case the measured quantities are the values 
of 

A,{z) = H{z)r,{za) (43) 



aX. z = 0.24 and z = 0.43 [59|, respectively. 



A^(0.24) = 0.0407 ± 0.0011 ± 0.0007, A^(0.43) = 0.0442 ± 0.0015 ± 0.0009 (44) 

with statistic errors and systematic errors. In this case, the parameter Qho is involved once 
again. 

Fitting results show that obviously this set of data has poor constraints on fi^o and gives 
poor improvement of FoM, as shown in Table II. 



III. COMBINATIONS OF DATA SETS 



In general, to get good constraints of DE parameters one often combines several types 
of observational data sets to fit cosmological models. In this section we study the effect of 
different combinations of the data sets introduced in the previous section on constraining 
DE models. We divide the combinations of data sets into two types: 
Type I: based on SN+BAOI+CMB, where the parameter fi^o is involved. 
Type II: based on SN+A+R, which only depends on the distance measurement. The calcu- 
lation is much simpler than the case of Type I. 

The basic combinations SN+BAOI+CMB (Bi) and SN+A+R (Bu) are two widely-used 
combinations in the literatures. The data sets to be combined with two basic combinations 
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data 








y2 • 

A, own 


FoM 


SN+BAOI+CMB 


n O7Q+0.017+0.033 
— u.u ±u — u . u^o 


1 07+0.14+0.27 
J- ' _n HQ— n 1 Q 


n 9Q+0.46+0.84 

— U . 1 u — i- . uu 


544.131 


27.721 


SN+BAOI+CMB+GRB 


n oQi +0.016+0.031 
u.zoi_n ni9— n 09"^ 


1 r,7+0.14+0.28 
-•-•^ ' _n OS— n 1 Q 

— u.uo — u. -Ly 


r. 9Q+0.45+0.82 
—0.78—1.74 


567.309 


27.401 


SN+BAOI+CMB+HUB 


n O7C+0.017+0.030 

— VJ.U J-U — \J.\J.^O 


1 nQ+0-13+0.25 
-'-•u«->_n OQ— n iQ 


n ^9+0.52+0.91 

— L/ . Urt — J. . t y 


554.034 


29.495 


SN+BAOI+CMB+RBAO 


n 97/1+0.016+0.029 
^•■^ ' ^—f) nnQ— n ni q 


-1 n9+0.12+0.24 
— u . uo — u . ± o 


n nf;+0-48+o.82 

—0.71—1.52 


544.721 


32.455 


SN+BAOI+CMB+CBMF 


n 97Q+0.019+0.034 

u.z ( y_ri noQ— n 099 

— u.uuy — u. 


-1 r,7+0.15+0.28 
J- .u ( _n nQ_n 9n 


n 00+0.51+0.90 

• —0.93— 1.76 


587.345 


26.319 


SN+BAOI+CMB+CBMF+RBAO 


n 97C+0.013+0.025 
^ ' °— 012— f) 023 


1 nQ+O.12+0.26 
•^"—0 08—0 17 


n nq+0.52+0.93 

U.UO_Q yQ_j^ gy 


587.990 


30.803 


SN+BAOI+CMB+HUB+RBAO 


n 97O+0.014+0.026 

U.Z(0_Q QQg_Q Q20 


-■ r,r|+0.11+0.22 
-'-•^'-'-0.09-0.18 


n na+0-50+0.86 

^•^°-0.60-1.40 


554.317 


33.736 


SN+BAOI+CMB+HUB+CBMF 


n OQQ+0-013+0.025 
'-'•^^^-0.007-0.017 


1 r,/l+0-15+0.29 
-■-•^^-0.10-0.21 


^•'^'-'-0.85-1. 71 


600.747 


36.651 


SN+BAOI+CMB+HUB+RBAO+CBMF 


n oQn+O-Oii+0.021 

"-'•^^"-'-0.008-0.017 


1 nn+0-16+0.28 

-•-•^^-0.08-0.19 


p, 71+0.55+1.03 
^•'-"--0.86-1.73 


602.267 


38.815 



TABLE III: The best-fitted values with 68% and 95% C.L. errors of ^mo, wq, wi, and figure of 
merit for Type I combinations in the CPL model. 



are GRB, RBAO, HUB and CBMF data. In what follows, we further study the properties 
of each data set based on basic combinations and their effects on fitting results; compare 
combinations of Type I and Type II; and find out the combination with the largest FoM. 

A. Properties of Data Sets with Two Basic Combinations 

1. GRB 

In Table III and IV, it is obvious that the results of SN+BAOI+CMB+GRB (with GRB) 
and SN+A+R+GRB are almost the same as those of SN+BAOI+CMB and SN+A+R 
(without GRB), respectively. This is in accord with the results in the previous section. 
We have also tested other combinations that contain the GRB data, and arrived the same 
conclusion. As a result, we conclude that at present, adding GRB data to other observational 
data docs not give better constraints on the DE parameters. 

2. RBAO 

This set of data is used only in Type I combinations. The inclusion of 
RBAO data always increases the value of the best-fitted Wq. In Table III, com- 
paring results of SN+BAOI+CMB, SN+BAOI+CMB+HUB, SN+BAOI+CMB+CBMF, 
SN+BAOI+CMB+HUB+CBMF with results of those with addition of RBAO data, we can 
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data 




Wo 






FoM 


SN+A+R 


n 977+0.014+0.029 


1 ni +0-14+0.27 
J- •vJi_n 1 1 _n 99 


r, r,Q+0.72+1.18 

— U.Uy_n Rfi_i 47 


542.642 


27.930 


SN+A+R+GRB 


n O77+0.014+0.030 
1 ' _n 014—0 027 


-1 n-i +0.13+0.27 
-'-•>J±_ou_o.23 


n r,c-+0.68+1.15 
^■^"-^—0 69—1 51 


565.873 


28.086 


SN+A+R+CBMF 


rv 977+0.015+0.029 
' -0.014-0.027 


r, QQ+0.12+0.26 
^•^^-0.14-0.25 


r, 9/1+0.84+1.35 
^•^^-0.63-1.48 


585.775 


27.181 


SN+A+R+HUB 


n 97/1+0.014+0.029 
"-"•^'^-0.013-0.025 


1 ni +0.14+0.27 
-'^•^-'--0.10-0.21 


n ni +0.56+0.99 
"-'•"-'-'■-0.73-1.53 


550.473 


30.316 


SN+A+R+HUB+CBMF 


n 97C+0.014+0.028 
"-"•^'^-0.013-0.025 


1 r,9+0.14+0.27 
-•-•^^-0.09-0.21 


r, m +0.53+1.00 
"-'•^^-0.78-1.58 


593.736 


29.519 



TABLE IV: The best-fitted values with 68% and 95% C.L. errors of firnO) '"^O) '"^i) and figure of 
merit for Type II combinations in the CPL model. 

see that adding RBAO data improves FoM well with only small increase of x^. Therefore it 
is worthwhile to contain RBAO data in the combinations to constrain DE models. 
3. HUB and CBMF 

From Table III and IV it is seen that containing HUB data can help to improve FoM, 
though the increase of FoM is much less than that of since there are 12 datapoints in the 
HUB data set. 

For CBMF data, there are 42 datapoints. With addition of this set of data, there is 
a large increase of x^, as seen in Table HI and IV, while the FoMs for combinations with 
CBMF are always smaller than those without CBMF, except the case that the HUB data 
is also included in Type I combinations. As a result, it seems not necessary to include the 
CBMF data to fit DE models. 

But it can be seen from Table HI that based on the basic combination SN+BAOI+CMB, 
the inclusion of both HUB and CBMF data can improve FoM well. For example, the 
FoM for SN+BAOI+CMB+HUB+CBMF is much larger than those for SN+BAOI+CMB, 
SN+BAOI+CMB+HUB and SN+BAOI+CMB+CBMF. While the data of HUB+CBMF 
have no similar effect in combinations of Type II, as shown in Table IV. Therefore, in 
Type I combinations it had better include both HUB and CBMF data, while for Type II 
combinations it is not necessary to contain the CBMF data. It is also shown that with data 
of HUB and CBMF the best fitted Q,mG is always larger than the ones in the cases without 
HUB or CBMF. 
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B. Comparing the Two Types of Combinations 

Both basic combinations of Type I and II consist of data from Type la Supernovae, 
Baryon Acoustic Oscillation, and Cosmic Microwave Background Radiations. The FoM and 
Qjno from two basic combinations (SN+BAOI+CMB and SN+A+R) are very close, and Xmin 
from SN+BAOI+CMB is a bit larger than that from SN+A+R since there is a bit more 
datapoints in SN+BAOI+CMB. In general, Xmin Type II are always smaller than those of 
the corresponding combinations in Type I. The fitting procedure of Type II combinations is 
much simpler than that of Type I. In Type II combinations all datapoints are only dependent 
on distance measurement (i.e., one needs only to use H{z)). In addition, the best-fitted Wq 
{^mo) in Type I combinations are always smaller (larger) than the corresponding ones in 
Type II combinations. 

C. Recommended Combinations 

Now we take the value of FOM to be the diagnostic to quantify the ability to constrain 
the DE parameters for different combinations of observational data sets. From Table III 
and IV, we can see the combinations with the largest FoM for two types of combinations are 
respectively: Combination I (COM I), SN+BAOI+CMB+HUB+RBAO+CBMF for Type 
I, with FoM=38.815 and Combination II (COM II), SN+A+R+HUB for Type II, with 
FoM=30.316. 

Though the FoM of COM II is much smaller than that of COM I, its fitting procedure 
is much simpler than that of COM I. The combination with the best FoM in Type I is just 
the one that contains the most data sets, while it is not the case in Type II combinations. 
The FoM for SN+A+R+HUB+CBMF (which owes the most sets in Type II combinations) 

is smaller than that for COM II (which does not contain the CBMF data). As shown in 
the previous section, adding CBMF data decreases FoM, and unlike the case in the Type 
I combinations the combination of CBMF and HUB data cannot improve FoM for Type II 
combinations well. 

To get tight constraints of DE parameters, we therefore suggest adopting COM I to fit DE 
models. On the other hand, if one prefers to use less data, a simper procedure and quicker 
calculation, COM II would be a good choice. Of course fitting results by using COM I or 
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COM II (or other combinations) will be a bit different. For the CPL model, the best-fitted 
Qmo from COM I is a bit larger than that from COM II. 

IV. SUMMARY 

We have analyzed the fitting properties of 11 widely-used observational data sets and 
their combinations. All data sets have been fitted with the CPL model, and the figure of 
merit (FoM) has been calculated for each case to estimate constraints on Wde- The FoM is 
considered as the diagnostic to quantify the abihty to constrain the DE model for different 
combinations of observational data sets. 

We first studied the constraint properties of 11 observational data sets and their combi- 
nations. The main results are: 

1. The case with the SN data used only gives a largish best-fitted flmo, and a relative 
small FoM, which means that the constraints on the DE parameters are very weak in this 
case. 

2. The addition of HUB data to SN data can help to increase FoM and the best-fitted 
Qmo becomes much smaller than that from the case with SN data used only. 

3. The addition of GRB data has almost no improvement in FoM, and fitting results of 
SN-I-GRB have almost no difference from those of SN. In this sense it seems not necessary 
to add GRB data in order to constrain simple DE models like the CPL model. 

4. The addition of 42 CBMF seems not quite useful to constrain DE models. 

5. With addition of CMB data to SN data, the FoM and constraints on the DE parameters 
are greatly improved than the case with SN data used only. 

We also compare fitting results with several BAO data sets. With these results, we 
studied different combinations of observational data sets and divided them in two types, 
based on two widely-used basic combinations: SN-I-BAOI-I-CMB and SN-I-A-I-R. We reach 
the following conclusions: 

1. It is helpful to constrain DE models by adding the RBAO data into combinations of 
Type I. 

2. In Type I combinations the inclusion of both HUB and CBMF data can improve FoM 
well, while in Type 11 combinations including CBMF data seems not helpful. 

3. The best-fitted Wq {flmo) in Type I combinations are always smaller (larger) than the 
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corresponding ones in Type II combinations. 

4. To get tight constraints of DE parameters, which means a larger FoM, we suggest to 
use the SN+BAOI+CMB+HUB+RBAO+CBMF (COM I) to fit DE models. On the other 
hand, if one prefers less data, simper procedure and quicker calculation, the combination 
SN+A+R+HUB (COM II) is a good choice. 

Obviously COM I and COM II will lead to a bit different fitting results of the cosmological 
parameters. For example, for the CPL model the best-fitted Qmo of COM I is larger than 
that of COM II. In some cases, different combinations of data sets give rather different fitting 
results. This is mainly due to the fact that at present some experiment data are not accurate 
enough and sometimes there may exist tensions among them. Obtaining tight constraints 
of Wde and finding the nature of DE need more accurate experiments. The method of FoM 
used in this paper can also be applied to more new observational data Q, Q H to find 
good combinations of data. 
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